Restless pions from orbifold boundary conditions: 
an explicit construction for noise reduction in lattice QCD 
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The exponentially decreasing signal to noise ratio in multibaryon correlators is the main obstacle 
to a first principles, QCD-based calculation of the nuclear force. Recently, we have proposed an 
orbifold boundary condition ("restless pions") that can dramatically improve this matter. Here we 
develop the idea further by proposing an explicit algorithm that can be used with purely periodic, 
"off the shelf gauge configurations. We also discuss finite volume corrections with the new boundary 
conditions and the use of the "Luscher formula" for the phase shifts. 



I. INTRODUCTION 



A first principles, QCD-based understanding of the origins of the nuclear force is still lacking. Besides the obvious 
theoretical interest, such a calculation would open the door for first principles calculations of other multi-baryon 
' quantities whose values are unknown and/or difficult to be measured and that impact nuclear, particle and astro- 
physics. A step forward in this direction was taken by recent unquenched calculations of baryon-baryon low energy 
fS i' phase shifts [JQ- They revealed what is currently the main difficulty facing this program: the exponential decrease 
D \ with (euclidean) time of the signal to noise ratio in multi-baryon correlators. We suggested recently 0] a method to 
, ■ alleviate this problem. It can be viewed as the use of a special kind of boundary condition along the spatial directions 
generated by an orbifold condition on the quark fields. One drawback of the method is that it leads to a non-positive 
fermion determinant. Thus, the only way it can be used in practice is by combining the new boundary conditions 
^ ■ on the valence quarks while keeping standard periodic boundary conditions for the sea quarks. While this raises 
, some conceptual issues, from the practical point of view the use of mixed boundary conditions is convenient as the 
' most costly part of a numerical lattice calculation - the generation of gauge configurations - does not have to be 
redone with the new boundary conditions. The purpose of the present paper is to make explicit the algorithm using 
previously generated periodic lattices and to discuss several issues raised, such as finite volume corrections and the 
validity of the Luscher formula [§, [l^, [HI ■ 

The origin of the noise increase with time in lattice calculations is well understood 0, 3 ■ A simple analysis shows 
that the signal-to-noise ratio of a two-nucleon correlator is given by 



where C{t) — {q^ {t)q^ {t)^^ {0)(f {Q)) is the two-nucleon correlator (spin, isospin and spatial indices are omitted for 
simplicity), cr^(t) the variance of C{t) after N measurements. Em and Et, are the lowest energies allowed for the 
nucleon and pion in the lattice. With periodic boundary conditions zero momentum states for both the nucleon and 
the pion exist and the minimum energy equals their rest mass, Em = M, E^r = rriTr- The basic idea employed in Q 
is to find a boundary condition where the pion zero- modes are forbidden resulting in Em = M but i?^ > to^. It is 
not obvious how to accomplish this since only the boundary conditions of the quark and giuon fields, not the hadron 
fields, are chosen at will. The solution found in 0] is to formally extend the lattice in one direction from < z < L 
to — L < z < L and impose the parity orbifold condition 
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q{x,y,-z,t) = V^q{x,y,z,t), 
A3{x,y,-z,t) = -A3{x,y,z,t), 

A^,{x,y,-z,t) = Af,{x,y,z,t), ^7^8, (2) 

where Vz = ij^'y^ implements a z-axis reflection. In other words, the fields in the z < region are not independent 
variables but merely the parity reflection of the fields in the 2: > region. The orbifold conditions for the quark and 
gluon fields imply in orbifold condition for the pion fields 



■^{x, y, -z, t) = -tt{x, y, z, t). (3) 

Equation ^ forbids the pion zero mode 7r(x, y, z, t) = constant. 

Since the fields in the z < are related to the ones in the z > region by a symmetry of QCD (parity), the action 
in both regions is the same, for every field configuration. Thus, there is no need to actually double the size of the 
lattice in an actual simulation. The only effect of the orbifold condition in equation ([2]) is right at the boundary at 
2 = 0. The precise form of the action at the boundary depends on the particular regularization being used. For 
instance, if we index the sites in the z-direction by z = 1, • • • ,L and z = —L, ■ ■ ■ , — 1 (no z ~ Q site!) and Wilson 
quarks are used we have 



S = K [9-1(73 - r)(7i - 91(73 + r)q_i] + (giqi + (7_iq_i) H 

= -2^91(73 + r)Vzqi + 2q^qi + • • • , (4) 

where k, is the hopping parameter, the index on the quark fields denotes the position in z (the remaining coordinates 
are implicit) and the dots denote the contributions from the two sides of the bulk, z > and z < (which arc equal 
to each other). We see then that the orbifolded [—L.L] lattice is equivalent to a [1, L] lattice with some extra terms 
residing at the boundary, as is the case with any boundary condition.^ Despite the fact that the z < region is not 
actually present in the simulations, we will use the extended lattice as a theoretical device in order to better analyze 
the properties of our construction. 

The orbifold condition of equation ^ (that we call "one-dimensional orbifold") is inconvenient because the spin 
up nucleon field zero modes are also eliminated. A slight change (that we call "three-dimensional orbifold") leads to 
a more manageable method: 



q{-x, 


-y, 


-z,t) 


= 'Pq{x,y,z,t), 


Ai{-x, 


-y, 


-z,t) 


= A4{x,y,z,t), 


A^i{—x, 


-y, 


-z,t) 


= A^{x,y,z,t), ^ ^ 0, 



where 7^ = 7^ implements a parity reflection. The pion and nucleon fields satisfy 



7r(-a;, -y, -z,t) -7r(x, y, z, t), 

N{-x,-y,-z,t) = N{x,y,z,t), (6) 

where in the last line we assumed the nucleons to be non-relativistic and used the non-relativistic limit 7^ ^ 1. 

As mentioned above, the orbifold boundary conditions cannot be used for the sea quarks due to the non-positivity 
of the determinant. Our proposal avoids this issue by using mixed boundary conditions: orbifold boundary conditions 
for the valence quarks and periodic for the sea quarks (and gluons). That means that periodic gauge configurations 
can (and have to) be used. That raises the question of how the orbifold condition can meld with periodicity. We will 
analyze this question in section II for both the one-dimensional and the three-dimensional orbifold and determine the 
allowed pion modes. In section III we will discuss the finite volume corrections in the orbifold construction and the 
use of the L'uscher formula. 



^ Notice that, contrary to the continuum case, the boundary conditions in lattice field theory are already contained in the action. Different 
lattice action terms localized at the boundary imply different boundary conditions in the continuum and the relation between them is, 
in general, a complicated dynamical question. 
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II. EXPLICIT LATTICE CONSTRUCTION 
A. "One-dimensional" parity orbifold 

For clarity, in this section we considered the simpler but less useful construction where the new boundary condition 
is imposed in only on direction. 

Consider a gauge configuration generated using standard periodic boundary condition for both gluons and (sea) 
quarks and index the sites on the z-direction by z = 1, 2, • • • , L. We can formally extend this configuration to one 
defined for z = f , 2, • • • , L and z = —L, • • • , — 1 by defining the gauge hnks in the z <0 region by 

Uz{x,y,z,t) = Ul{x,y,-z,t), 

U^,{x,y,z,t) = U^,{x,y,-z,t), for pi ^ 3. (7) 

An ensemble of gauge configurations generated by this procedure is identical to one that would be generated by 
actually simulating a larger —L<z<L lattice where the fields in the two halves are constrained to obey the orbifold 
condition: 

qs{x,y,-z,t) = Vzqsix,y,z,t), 
A3{x,y,-z,t) = ~A3{x,y,z,t), 

A^,{x,y,-z,t) = Af,{x,y,z,t), ^7^3, (8) 

where qs are the (sea) quark fields. Notice that there is no gauge link connecting the z = 1 and z = —1 sites (see 
figure [TJ 




FIG. 1: Lattice where the gauge and sea quarks "live". Only the z-coordinate is shown. The fields on the left sub-lattice are 
not independent but only a parity reflection of the ones on the right sub-lattice. 

The action for the valence quarks is also defined on a doubled lattice. The difference from the gauge and sea quark 
sector is that hopping terms arc added connecting the z = 1 to z = —1 and z = L and z = —L for all values of x, y 
and t. In these links the gauge field is defined to be trivial {Uz = 1). This construction is pictured on figure [2l The 
valence quark fields at z < are, again, not independent of their counterparts at z > 0; qy{x,y, —z,t) is constrained 
to be the parity reversal of qv{x, y, z, t) 

qv{x,y,z,t) ::^'Pzqv{x,y,z,t). (9) 

Notice that, in the valence action, the nodes with z ~ 1, z ~ L have one more nearest neighbor that the remaining 
sites. 

The new plaquettes in the action involving the links between z = I and z = — I (z = L and z = —L ) arc trivial 
and don't contribute to the action. The action in the z < sub-lattice, for any field configuration, is the same as 
the one in the z > sub-lattice, on account of the invariance of the lattice QCD action under parity transformations. 
Thus, in an actual simulation, only the z > sub-lattice is used. The only effect of our construction in the valence 
action is the addition of an extra term at the boundary at z — I. The precise form of this term is determined by the 
particular lattice action used and the orbifold condition in eq. For the Wilson action, for instance, we have the 
result in equation 

In the continuum limit, the orbifold conditions in cq. (|8l9p imply in 
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FIG. 2: Lattice where the gauge and valence quarks "hve". Only the z-coordinate is shown. The fields on the left sub-lattice 
are not independent but only a parity reflection of the ones on the right sub-lattice. The double line denotes a hopping term 
with trivial gauge field f/z = 1- 



qs{x,y,-z,t) = Vzqs{x,y,z,t), 

qy{x,y,-z,t) = Vzqv{x,y,z,t), 

Af^{x,y,-z,t) = Af,{x,y,z,t), for 3 

A3{x,y,-z,t) = -A3{x,y,z,t), ior (10) 

The sea quark fields are periodic along the z-direction (going from z — +0 to z = L or from z = —0 to ^ = —L) but 
can be discontinuous in going from z = +0 to z = —0 since there is no hopping term connecting the two halves of the 
extended lattice. The valence quarks, however, are also continuous along the z — +0 to 2; = —0 border. 

We can now use equation (|10p to figure out the boundary conditions satisfied by the hadrons in the continuum 
limit. In the spirit of partially quenched chiral perturbation theory @, 0, there are three kinds of pions, depending 
whether their constituents quarks and antiquarks are sea or valenve quarks. For all three kinds of pions we have the 
parity orbifold condition: 

iTyyix^y^ z^t^ '^vvi^x J y ^ z ^ ^ 

Trysix,y,-z,t) = -Hys{x,y,z,t), 

■^ss{x,y,-z,t) = -T^ss{x,y,z,t). (11) 

Furthermore, the translation invariance in the sea sector implies the periodicity of the ttss field in all three spatial 
directions. That is, not only iTssi+O^y, z,t) = i:ss{L,y, z,t), Trss{—0,y,z,t) = i:ss{—L,y,z,t) (and similarly in the 
y and z directions) but its derivatives are also the same ( 9a;7rss(+0, y, z, t) = dxT^ss{L,y, z,t), ■■■ ). However, 
7rss(-|-0, y, z, t) is not necessarily the same as 7rss(— 0, y, z, t). Thus, the condition iTssix, y, —z, t) = —iTssix, y, z, t) does 
not imply in tt^s {x,y,0,t) = 0. 

On the other hand, the continuity of at z = (the result of the extra links) implies in 

'Kvv{x,y,L,t) = 'Km,{x,y,+Q,t) = ~'Kyy{x,y,~Q,t) 'Kyy{x,y,-L,t) 0. (12) 

The vanishing of 7r„i, at z = eliminates its zero mode which is the main point of our construction. ■Ugy satisfies the 
same boundary conditions as Uss- 

The boundary conditions discussed above specify the allowed modes for each kind of pion: 

^™(z) ~ e''-r(nMyy) 71,, 71^ = 0, ±1, ±2, • • • , = 1,2,... (13) 

^.4z),^„4z) ~ e*T(".-+»«J.)|^^g^W^^^^^ , n„n, = 0,±l,±2,..., = 0, 1, 2, • • • . (14) 

The discontinuity for tt^s, t^sv may seem strange. Remember, however, that there is no link between z > and z < 
in the sea quark action, so the discontinuity costs no energy. It is important to note that the presence of zero modes 
for TTys and TTss does not spoil the noise reduction of baryon correlators. This is because the sea-sea and sea-valence 
pions can only appear in intermediate states addition to the valence pions. ^ 



^ Notice that in correlators of two (valence) nucleons will always include, at least, 6 valence pions and thus the minimal energy of allowed 
intermediate states will be bounded from below by (six times) the minimal energy of the valence pions. 
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The valence nucleons (those composed of three valence quarks) are periodic in the x and y directions, continuous 
at z = and satisfy the orbifold condition N{x, y, — z, t) = VzN(x, y, z,t). At z = and in the non-relativistic limit 
this implies in 

Nix,y,0,t)^<j,Nix,y,0,t). (15) 
The up and down spin components obey different boundary conditions. 



B. "Three-dimensional" parity orbifold 



This case is very similar to the "one-dimensional" orbifold case. Again, we start from gauge configurations periodic 
in all three directions. This time, however, we choose the x and y coordinates to vary from —L/2 to L/2 while the z 
coordinate varies from 1 to L. We formally extend the gauge configurations to the z < region through the relations: 



qs{~V,t) - Vqs{r,t), 
A{-r,t) = -A(r,t), 

Ai(-r,t) = A4(r,t), (16) 

with T' = 74 is the usual parity operator. The valence quark propagators are computed in this background after 
trivial {U = 1) links are added connecting {x,y,l,t) to (— .t, — y, — 1, t). Again, the result of having the extended 
lattice is simply an extra term in the lattice action living at the first node z = 1 whose precise form depends on the 
fermion discretization used. Contrary to the "one-dimensional" orbifold case though, this extra term is non-local and 
corresponds to hops from {x,y, l,t) to/from {—x, —y, l,t). This construction is equivalent as having valence quark 
fields that are periodic in the a;, y and z directions and obey the orbifold condition 

g,„(-r,t) =7'<z„(r,t). (17) 

In addition (and contrary to the sea quarks), is continuous in going from (x, y, +0, t) to (— x, — y, —0, t). Using the 
Wilson lattice action as an example, an explicit construction of the valence quark action, we have'^ 

S = K [qx.y,l,th3 + r)qx,y-l.t - qx,y-l,t{l3 " f )fe,y,l,t] + {qx.,y,l,t fe,y,l,t + qx,y~l.t -l,t) H (18) 

where as in eq. (|3|), the ellipses denote the action from the two sides of the bulk. Using the three-dimensional 
orbifolding conditions, eq. (jl7p . one can show that the action is given by 

S = K [qx,y,i,th3 + 'r)'Pq-x,-yS.t - q-x-yS.tVijs ~ r)qx^ys,t] ~ [qx,ys,t qx,y.i,t + q-x-y,ld q-x-y,i,t\ H 

= 2k qx,y,l.t{l3 + 'r)'Pq-x,-y,l,t - '2qx.,y.l,t qx,y,l.t H (19) 

where the second line comes from a relabeling of the dummy spatial indices under summation. Notice the dif- 
ferences between this boundary action and that for the one-dimensional orbifolding construction, eq. ([4]). In the 
one-dimensional case, the extra boundary action is entirely local while in the three-dimensional case, eq. p9p . the 
boundary action is explicitly non-local. 

We can now determine the boundary conditions satisfied by the hadronic fields in the continuum limit. Like in the 
"one-dimensional" orbifold 



7r™(-r,t) = -7ri,t,(r,t), 
nys{-r,t) = -7r„s(r, t), 
nss{-r,t) = ~Trss{r,t). 

(20) 

The valence pions. besides obeying the same boundary/orbifold conditions as the sea and sea-valcncc pions, are also 
continuous along the z = plane. This implies on the vanishing of ir^y at the origin (0, 0,0,t) since 7r.u„(0, 0, 0, t) = 



^ Recall that there is no z = site, for simplicity of writing down the action. Also, there is a similar contribution to the action at z = L. 
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— 7r^„(0, 0, 0, t) = 0. Due to the periodicity, 7r„„ then vanishes at the border of the lattice as can be seen through the 
foUowing relations, 

^™(i/2,0,0,t) = -7r,„(-i/2,0,0,t) = -7r„„(L/2,0,0,0 = 0, (21) 

and similarly for the y and z directions. The first relation arises from the parity-orbifold construction while the second 
relation arises from the periodicity of the lattice. 

Taking into account the boundary conditions satisfied by the tt^v , i^sv and tt^s fields we find that the allowed modes 
are 

2nxT^x 2n„7ry n^Trz 
TTyy(x,y,z) ^ sm( — - — )sm( — - — )sm(— 



TTss{x,y,z),'!Tsv{x,y,z) ~ sm( — - — )sm( — ^ — ^ 



sign(z) cos(2^) 



rsign(z)cos(2^) 

eos(^)cos(^) 

isin(^) 

. 2n^7ra; 2nyiiy .n^Tiz 
sin{—j—) cos( ^ ) cos(— ^), 

cos(^^)sin(^^)cos(^). (22) 



Notice that the minimal energy of a valence pion is y^2(27r/L)^ + [tt/LY + rri^, a value larger than the one quoted 
in where the appropriate boundary conditions at z = L were not taken into account; see figure [3] for an updated 
estimate of the improvement in the signal to noise. 

Nuclcon fields satisfy, as a consequence of equation p?|) . the boundary condition N{—r,t) — 7^A^(r, i). In the 
non-relativistic limit this reduces to N{~r,t) — N{r,t). Thus, unlike pions, the nucleon zero mode is allowed. This 
feature is very convenient to construct two-nuclcon states at rest and use the Liischer formula to extract the zero 
energy scattering phase shift. 



2M-3 
1000 




FIG. 3: Estimated energy splitting which controls the decay of the signal-to-noise for the two nucleon system, 2M — SE-n. The 
(red) dashed line is for periodic boundary conditions with Ett = m-n while the blue line is for the three dimensional orbifold, 
with E-n = i/9(7r/L)2 -f m^. , in which we have held m^L = 4 fixed in our estimation. The pion mass dependence of the nucleon 
mass was estimated from Refs. [13, [3 • 



III. FINITE VOLUME CORRECTIONS 



The finite volume corrections to quantities computed in the lattice are, as a long distance property, dependent 
on the boundary conditions used. In periodic lattices, masses, decay constants and scattering lengths receive only 
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exponentially suppressed finite volume corrections proportional to ~ ^-m^L^ Unfortunately, the absence of pion zero 
modes enhances these finite volume corrections and they become only power law suppressed. It is simple to estimate 
these eff'ects. For observables describable by a low energy effective theory like chiral perturbation theory these effects 
are dominated by pion loops. The difference between a one-loop diagram computed at finite volume and in the infinite 
volume limit is of the form 

^^ = jiJ:m~J^m' (23) 

for some function f{q), where the sum is over the allowed momenta for the pions. In periodic lattices q = 21:71/ L with 
ft having integers components. In this case it can be shown that, as long as f{q) has no singularity along the real axis 
(that is, the particles in the loop cannot be on-shell simultaneously), ~ e~™"^. When the allowed momenta do 
not include the zero mode can be larger. Let us take, for instance, the case of the "three-dimensional" orbifold. 
The sum over the allowed valence pion modes is 



1 ^ , f^nn^ 2Tmy irnA f (Pq 



L3 ^ ■' \ L ' L ' L J J (27r)3 

- ^_ \^ f i 2^"^ ^TTTT,^ 7rn^ \ ^ [ d?q 

~2L3 ■' [ L ' L ' L J J {2^f^^'^' 

1 „/27rn^ 27rny ■nn^\ f d^q 1 

n^,ny,n^=0,±l,±2,--- ^ / J \ ) 

= O(e-™"^)-^/(0,0,0), (24) 

where we assumed that /(i?) is an even function. The /(0,0,0)/L3 term is the power law finite volume correction 
mentioned above. 

Through a more detailed analysis one can estimate the numerical value of these power law corrections. Since this 
analysis is dependent on the particular observable of interest we will refrain from pursuing it here. However, we note 
that because of chiral symmetry, the pions, at leading order, are derivatively coupled to hadrons and so these power 
law corrections will typically be further suppressed in the chiral power counting. We would like to further comment 
regarding the impact on the extraction of phase shifts by the use of the Liischcr method [1, |Tl[ , since the main 
application we have in mind for our construction is the study of nuclear forces. In the Liischer method, the energy 
levels of a two-particle system are related to the (infinite volume) phase shift. In the case of boxes sizes L much larger 
than the scattering length a between the particles this relation is 



. „ 47ra 
Ml\E ^ 



l + ci| + C2(|)2 + ... , (25) 



where C\ and C2 known numerical constants, Ai? = E — 2M the shift in energy between the non-interacting and the 
interacting two-particle state and M the mass of each particle. With spatially periodic boundary conditions, this 
formula is exact up to corrections suppressed by ~ e-^xi- xhe ellipsis in equation includes corrections 

due to the mixing of partial waves, starting at order In the case of the orbifold boundary conditions, the first 

two terms in equation (j25p are unchanged since the power law corrections to the energy levels described by it are 
due to the on-shell propagation of nucleons whose boundary conditions are the same as in the periodic case. The 
large (~ ^/L^) finite volume corrections to the mass of each particle cancel against the same correction to the energy 
levels of the two-particle state and does not affect the extraction of the scattering length. In other words, the finite 
volume corrections to E are cancelled by those of 2M [l^. There will be, however, corrections to the nucleon-nuclcon 
interaction (and consequently to the value of a) that are suppressed by l/L^ only. However, combining it with the 
l/L^ factor present in equation (j25|l . the finite volume correction to the energy level will appear only at order 
This is paranietrically smaller than several other corrections coming from contamination from higher partial waves, 
etc.. and arc not a cause for concern, at least as long the condition a ^ L is valid. The breaking of hypcrcubic 
invariance (due to the special role of the z direction) is more cause for concern as it allows for a mixing of higher 
partial waves beginning at order 
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IV. CONCLUSIONS 



We spelled out in detail an algorithm to implement orbifold boundary conditions in lattice QCD that eliminates the 
pion zero mode and, consequently, improves the signal-to-noise ratio in multi-baryon correlators. Our construction 
involves the use of a sea sector with periodic boundary conditions (and the use of standard, "off-the-shelf gauge 
configurations) and a valence sector with orbifold boundary conditions. In the most favored incarnation of our 
construction (the "three-dimensional orbifold" ) the energy of the lowest pion mode is given by •\/9(7r/L)2 -I- what 
is enough to completely eliminate the exponential growth of error in two-nucleon correlators for pion masses larger 
than about « 300 MeV and lattice sizes L « 4/TO7r (and drastically reduce it for hghter pions). We pointed out 
that the absence of the pion zero mode enhances finite volume corrections, such that they are of order ~ 1/^'^, as 
opposed to e"™"^ as in the periodic case. These enhanced finite volume effects, however, do not preclude the use of 
the Lilscher formula to extract baryon-baryon phase shifts. Indeed, the corrections to the scattering length extracted 
through the Liischcr formula are higher order in powers of 1/i and are furthermore suppressed in the chiral expansion. 
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